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Here we show how an anisotropic fluid in the diffusion limit can be equivalent to an isotropic 
fluid in the streaming out limit, in spherical symmetry. For a particular equation of state this 
equivalence is total, from one fluid we can obtain the other and vice versa. A numerical master 
model is presented, based on a generic equation of state, in which only quantitative differences are 
displayed between both non-adiabatic fluids. From a deeper view, other difference between fluids 
is shown as an asymmetry that can be overcome if we consider the appropriate initial-boundary 
conditions. Equivalence in this context can be considered as a first order method of approximation 
to study dissipative fluids. 
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I. INTRODUCTION 

Emission of massless particles is an important process 
when studying evolution of massive stars. If the fluid is 
nonadiabatic two transport mechanisms are recognized 
as extreme and opposite, namely, the streaming out and 
the diffusion limits. In the event of gravitational collapse 
these two transport regimes alternate, playing different 
roles depending on the configuration and physical pro- 
cesses [J, g], [3j. 

It is well known that different energy-momentum ten- 
sors can lead to the same spacetime [l]-[9]. Different 
kinds of physical phenomena may take place to give rise 
to local anisotropy [10 . The influence of anisotropy in 
self-gravitating systems has been studied mostly under 
static conditions (see [11] and references therein). In 
|12j a general study of spherically symmetric dissipative 
anisotropic fluids was reported, with emphasis on the re- 
lationship between the Weyl tensor, the shear tensor, the 
anisotropy of the pressure and the density inhomogene- 
ity. For instance, shear viscosity can be considered as a 
special case of anisotropy [T3]— [TS] , [TB], with no contri- 
bution to energy density and energy flux as heat flow. 
Bulk viscosity does not change the degree of anisotropy. 
Electrically charged fluids also can be seen as anisotropic 
fluids with no contribution to heat flow [IT], [15] . 

Many years ago the equivalence between radiating per- 
fect fluids and anisotropic fluids was shown |19j . which 
also applies under certain conditions to multicompo- 
ncnt fluids [20], [21] ■ Recently it was proposed that an 
anisotropic fluid in the diffusive limit is, in some sense, 
the most fundamental model [55] because it can absorb 
the addition of shear viscosity, electric charge and null 
fluids. In fact, the free streaming process generates heat 
flow, radial pressure (nonisotropically) and energy den- 
sity. In consequence the resulting fluid is anisotropic. 
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We show here how in practice an isotropic fluid in the 
streaming out limit can be equivalent to an anisotropic 
fluid under a diffusion transport mechanism. We argue 
that an isotropic fluid with free streaming is as funda- 
mental as an anisotropic fluid with heat flow. We il- 
lustrate our claim numerically by means of a master 
model. To do that we use the seminumeric method 
known as the postquasistatic approach (see [53] and ref- 
erences therein), which fits well for this paper purposes. 
We do not know of any other hydro solver in numerical 
relativity that considers a radiating fluid in either regime, 
streaming out or diffusion, isotropic or anisotropic, in 
spherical symmetry. 

The equivalence between nonadiabatic fluids as pre- 
sented here can be used to study different situations in 
the gravitational collapse with the same fluid, but this is 
not always possible. The equation of state seems to be 
crucial as well as the initial conditions, the exterior space- 
time and the junction conditions. The equivalence has 
been used as a method to solve the field equations, but 
not to extract physics. This situation is quite similar to 
a massless scalar field viewed as a fluid. The asymptotic 
conditions near the center of symmetry, at the boundary 
surface and near infinity are different for the scalar field, 
in comparison with a bounded source of matter [24] , [25j . 

The remaining parts of this work are organized as fol- 
lows: In Section II, we show how a local comoving ob- 
server can handle the equivalence between non-adiabatic 
fluids and how the field equations can be recast. In Sec- 
tion III, we give some examples using a master numerical 
model. In section IV, we revisit the results in the light 
of previous ones and conclude with some remarks about 
the limits and scope of the equivalence of nonadiabatic 
fluids. 



II. EQUIVALENT RADIATING FLUIDS 

We can clearly show at once the aforementioned equiv- 
alence between nonadiabatic fluids using Bondian ob- 
servers (comoving in a locally Minkowskian spacetime) 
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to describe how physics comes out from hydrodynamics. 
The covariant energy-momentum tensor for a radiating 
perfect fluid as seen by a comoving observer with respect 
to a local one with radial velocity U) is 




(1) 



representing a radiating isotropic fluid of energy density 
p, pressure p, and energy density e traveling in the radial 
direction. If we define 
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= p + e, 
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Pr 
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Pt 
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pressure pt and heat flow in the radial direction q. Ob- 
serve that the degree of anisotropy is 



A=pt-p r = -e, (4) 



which in turn is seen as an equation of state. As a mat- 
ter of fact, any comoving observer sees such equivalence, 
it does not depend on the system of coordinates. No 
rotation about the timelike two-plane is required [19j . 
[2"U] . [2"o] because they are implicit in and It is 
manifest that the covariant energy-momentum tensor is 
invariant. 

We have to write the field equations in one limit to 
show explicitly and formally how we can get the equiv- 
alence simply recasting them. Using Bondi's metric in 
spherical form [27J 



the fluid can be seen as 
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ds 2 = e 2/3 ( - du 2 + 2dudr) - r 2 (d9 2 + sin 6 2 dej) 2 ) , (5) 
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which can be interpreted as an anisotropic fluid with ef- 
fectives energy density p, radial pressure p r , tangential 



where (3 = /3(u, r) and V = V(u,r). The field equations 
for an isotropic fluid in the free streaming approximation 
can be written as [28] 
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where a comma denotes partial differentiation with re- 
spect to any coordinate, and the mass aspect m is defined 
by means of 
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As expected, we can write the LHS of the field equations 
in the following way, only rearranging terms 
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where p and p are the so-called effective variables. In 
standard numerical relativity S and p are the conserva- 
tive variables, and p the flux variable. Now (|6])-(|9| are 
the field equations for an anisotropic fluid in the diffu- 



sion limit with the LHS formed with (!!)-( 13 1, &2dh and 
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the additional equation of state Q. In this sense it is 
said that an isotropic fluid in the free-streaming limit is 
equivalent to an anisotropic fluid in the diffusion limit. If 
the fluid is described ah initio as anisotropic in the diffu- 
sion limit and the equation of state to take into account 
the degree of anisotropy is given by Q , then we get nec- 
essarily an isotropic fluid with free streaming. We shall 
explore a generic situation, that is, with any other equa- 
tion of state. In that case: does the equivalence between 
these two nonadiabatic fluids hold? 

In order to gain physical insight we write the field equa- 
tion |9]) as the generalized equation for hydrostatic sup- 
port of Tolman-Oppenheimer-Volkoff [35] , [35] 
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One way to illustrate our view point is using a method 
which represents the departure from equilibrium in the 
gravitational collapse problem [28 . That method re- 
ported many years ago could be used as a testbed for 
numerical relativity |30) . It as been interpreted as a 
postquasistatic approximation [23]. The procedure is 
based on the cornerstone paper of Bondi to deal with 
radiating matter [27] . 

For details about the method used to solve the equa- 
tions, matching conditions and its connection with stan- 
dard numerical relativity see [25J , [23] , [30] ■ Here we only 
write the appropriate boundary condition |31j 



[Pr]a = qa, 



(15) 



which is a consequence of the Darmois-Lichnerowicz con- 
ditions, that is, the continuity of the first and second 
differential forms: 



m = M(u), 



/3a =0, 
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where the subscript a indicates that the quantity is evalu- 
ated at the surface r — a(u), rh a = m and M corresponds 
to exterior Vaidya spacetime 




FIG. 1: Evolution of the radius a for different models: Ml; 
M2; M3- with h = 0.33; M3+ with h = 1.33; M4- with h = 
0.33; M4+ with h — 1.33; all radiating the same fraction 
of the initial total mass (1%). The initial conditions are: 
A(0) = 5; F(0) = 0.6; £1(0) = 0.8333. Model M4+ may 
replicate model Ml if the anisotropy parameter h w 1.2. 
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ds 2 = f 1 - — ) du 2 + 2dudr - r 2 (dd 2 + sm9 2 d(j) 2 ). 



Clearly in the case of free streaming the boundary con- 
dition reads 



FIG. 2: Evolution of the energy density p (multiplied by 10 3 ) 
for the streaming out limit model (Ml) and initial conditions 
A(0) = 5; F(0) = 0.6; £1(0) = 0.8333. 
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To numerically study the equivalence, we shall use a 
Schwarzschild-like model with intrinsic anisotropy based 
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on the following generic-equation-of-state (GEoS): 



F=l- 



2M 



(r — 2m) 



(21) 



where C is a constant. This GEoS has been used in the 
past for a number of physical situations which represent 
phase transition, viscous and electrically charged fluids, 

E2, na, Dai, \m, ea- 



rn. MASTER MODEL 

Here we describe a master model, which includes all the 
possible configurations of interest for our purposes, based 
on the models described in previous investigations [2"5] . 
[31], [29], [13]. Following the postquasistatic protocol, 
we build a Schwarzschild-like master model which cor- 
responds to an incompressible anisotropic/isotropic fluid 
as the static "seed" interior solution. Thus, the effective 
energy density and the effective pressure are 
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With these effective variables and integrating ^ and Q 
we obtain 
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The system of equations at the surface is 
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overdot denotes the total derivative with respect to r, L 
is the luminosity as a function of time given in some con- 
venient way. Thus, the master model summarizes in one 
equation, (26c), four previously reported models of the 



same type (Schwarzschild-like) . The whole model is rep- 
resented by the set of equations (26) and allows to control 
the transport mechanism and anisotropy using only two 
parameters (h and I), h = 1 represents isotropy, h < 1 
implies pt > p r , and h > 1 means p t < p r . On the other 
hand, I = 1 is for the streaming out limit and i = for 
the diffusion limit. Therefore, this master model contains 
all possible situations, that is, isotropy in the streaming 
out limit (Ml: h = 1; £ = 1); isotropy in the diffusion 
limit (M2: h = 1; £ = 0); anisotropy in the streaming 
out limit (M3±: h ^ 1; I = 1); anisotropy in the diffu- 
sion limit (M4±: h ^ 1; £ = 0). Equivalent models come 
from the corresponding master models depending on the 
equation of state used, for instance, El comes from Ml 
usin g p] ); E2 comes from M4+ which in turn is built us- 
ing ( |21[ ). In other words, equivalent models are obtained 
with no additional calculations but just singling out the 
source model data. 

With the following set of initial conditions we can inte- 
grate numerically the system of equations at the surface 
and go inside the distribution without restrictions: 

A(0) = 5; F(0) = 0.6; fl(0) = 0.8333 

with L given as a Gaussian pulse carrying away a fraction 
of total initial mass (w 1%): 

L = i e-(---o)7< 

where a = A and To = 5. Figure 1 displays a typical situ- 
ation in which we can see how the generic fluid behaves 
depending on anisotropy and transport mechanism. Al- 
though the behavior shown in figure 1 for each model has 
been reported in previous investigations [28], |31j . [29) . 
[13] , for the sake of completeness and to gain physical 
intuition we summarize and explain these results here. 
It is clear from Eq. ( 14 ) that anisotropy, second term 



m(0)' 



of the RHS, increases (or decreases) the effective grav- 
itation depending on the tangential pressure and irre- 
spective of the transport mechanism [29 . On the other 
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FIG. 3: Evolution of the radial pressure p (multiplied by 10 4 ) 
for the streaming out limit model (Ml) and initial conditions 
A(0) = 5; F(0) = 0.6; fi(0) = 0.8333. 



+J 

■rH 

U 
O 

iH 

> 



(0 

o 
o 



(0 

■rH 

TS 
(0 
PS 



0.05 



0.00 



-0.05 



-0.10 



-0.15 



-0.20 



-0.25 



r/a=0.2 
r/a=0.4 
r/a=0.6 
r/a=0.8 
r/a=1.0 



2 4 6 



8 10 12 14 16 18 20 
Time 



FIG. 5: Evolution of the radial local velocity uj for the 
streaming out limit model (Ml) and initial conditions A(0) = 
5; F(0) = 0.6; «(0) = 0.8333. 
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FIG. 4: Evolution of the energy flux e (multiplied by 10 3 ) 
for the streaming out limit model (Ml) and initial conditions 
A(0) = 5; F(0) = 0.6; Q(0) = 0.8333. 



FIG. 6: Evolution of the energy density p = p + e (multiplied 
by 10 3 ) for the equivalent diffusion limit model (El) and ini- 
tial conditions A(0) = 5; F(0) = 0.6; n(0) = 0.8333. 



hand, a positive heat flow always diminishes the effective 
gravitation irrespective of the anisotropy [31] • For these 
reasons the radius for model M4+ (p t < p r ) is smaller 
than that of model M4— (p t > p r ); model Ml (p t = p r ) 
is somewhere in between. Figures 2-5 display results for 
model Ml and figures 6-7 for the equivalent model El. 
Note that we get p and p r of model El from p, p and e 



of model Ml; all the other physical variables (q, lu and 
p t = p) are undistinguishable for these two equivalent 
models. Finally, making numerical experiments we find 
that model M4+ behaves at the surface like model Ml 
when h — 1.2, which is shown in figures 8-12. The equiv- 
alent model E2 is displayed in figures 13 and 14. Observe 
how physical variables in model El (figures 6 and 7) be- 
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FIG. 7: Evolution of the radial pressure p r = p + e (multi- 
plied by 10 4 ) for the equivalent diffusion limit model (El) and 
initial conditions A{0) = 5; F(0) = 0.6; Q(0) = 0.8333. 



FIG. 9: Evolution of the radial pressure p r (multiplied by 
10 4 ) for the diffusion limit model (M4+) with h = 1.2 and 
initial conditions A(0) = 5; F(0) = 0.6; fi(0) = 0.8333. 
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FIG. 8: Evolution of the energy density p (multiplied by 10 ) 
for the diffusion limit model (M4+) with h = 1.2 and initial 
conditions A{0) = 5; F(0) = 0.6; Q(0) = 0.8333. 
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FIG. 10: Evolution of the heat flow q (multiplied by 10 3 ) 
for the diffusion limit model (M4+) with h = 1.2 and initial 
conditions A{0) = 5; F(0) = 0.6; 0(0) = 0.8333. 



have like physical variables in model M4+ (figures 8 and 
9) and how physical variables in model E2 (figures 13 
and 14) behaves like model Ml (figures 2 and 3). Ex- 
cept for the tangential pressure, the isotropic fluid in the 
free streaming limit can be equivalent to the anisotropic 
fluid in the diffusion limit. The GEoS works remarkably 
well to illustrate the equivalence, displaying only quan- 



titative differences for models El and M4+; E2 and Ml. 
In this sense, the two equivalent models El and E2 are 
qualitatively the same as M4+ and Ml, respectively. It 
suggests that if we impose \p t ] a = the performance of 
the GEoS will get better, but this extra condition is out 
of our master model. 
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FIG. 11: Evolution of the local radial velocity cu for the dif- 
fusion limit model (M4+) with h = 1.2 and initial conditions 
A(0) = 5; F(0) = 0.6; fi(0) = 0.8333. 



FIG. 13: Evolution of the energy density p — p — q (multiplied 
by 10 3 ) for the diffusion limit model (E2) with h — 1.2 and 
initial conditions A(0) = 5; F(0) = 0.6; 0(0) = 0.8333. 
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FIG. 12: Evolution of the tangential pressure p t (multiplied FIG. 14: Evolution of the radial pressure p — p r — q (multi- 
by 10 4 ) for the diffusion limit model (M4+) with h = 1.2 and plied by 10 4 ) for the diffusion limit model (E2) with h = 1.2 



initial conditions A(0) = 5; F{0) = 0.6; Q(0) = 0.8333. 



and initial conditions A(0) = 5; F(0) = 0.6; fi(0) = 0.8333. 



We want to close this section with the following com- 
ment. Setting up the system we realized that for a com- 
pact configuration, initially at rest, the fraction of radi- 
ated mass was bounded in the diffusion approximation. 
Otherwise the local radial velocity was out of the real 
domain. This situation changes radically with a greater 



initial velocity of collapse or a lesser fraction of radiated 
mass (less than 0.01%). For the free streaming this re- 
striction does not exist, allowing to radiate of more than 
a 10% of the initial total mass. This "asymmetry" is not 
in favor of general equivalence but this was obviously 
overcome. 



IV. CONCLUDING REMARKS 

We have shown how radiating fluids can be equiva- 
lent, when the interior transport mechanisms are oppo- 
site and extreme, namely, for the free streaming and the 
diffusion approximations. From the formal point of view 
the equivalence between nonadiabatic fluids is indistin- 
guishable when applied to an isotropic fluid with free 
streaming and to a diffusive anisotropic fluid. One can 
be derived from the other. In order to show numerically 
this equivalence, without loss of generality, we have used 
the postquasistatic approach to consider fluids leaving 
equilibrium. The equivalence can be total or partial, de- 
pending on the initial-boundary conditions and the equa- 
tion of state. We illustrate by means of a master model 
how an isotropic fluid with free streaming can be equiv- 
alent to an anisotropic with heat flow dissipation. They 
are mathematically and physically equivalent from each 
other when using the equation of state Q. In a generic 
situation the equivalence depends strongly on the equa- 
tion of state, initial conditions, junction conditions and 
exterior spacetime. We have used in the past the GEoS 



(21 1 to model different physical scenarios. Viscous [16] or 
electrically charged |T5] fluids can be seen as anisotropic, 
but they are not equivalent to isotropic fluids with free 
streaming or anisotropic fluid with heat flow. 

When coherent radiation is important (scalar, electro- 
magnetic or gravitational radiations) the associated flu- 
ids induce anisotropy and heat flow but the observers are 
always resting at infinity. 

Within the family of fluids than can be equivalent it 
is then plausible that what anisotropy gives (stiffening) 
heat flow takes (softening) . In fact, in this case the ex- 
terior spacetime is the same from the geometrical point 
of view. But a judicious observer realizes that, in gen- 
eral, dissipation in the streaming out limit can never be 
as that of heat flow. Mathematically, the reason can be 
easily understood. In the streaming out limit the alge- 
braic equation to get the radial local velocity from the 
field equations is a linear expression, while in the diffu- 



sion limit the analogous equation is quadratic. In con- 
sequence the associated discriminant imposes a severe 
restriction on the initial velocity and the fraction of total 
mass that can be radiated away. We have shown that the 
only exception emerges when the two equivalent fluids are 
linked by the particular equation of state Q . Apparently 
the diffusion mechanism does not accept static equilib- 
rium. A complex velocity is as unacceptable as a negative 
energy density or a superluminous local radial velocity. 
Therefore a complex velocity is not physically reliable. 
Interestingly, in comoving coordinates that asymmetry 
is not present anymore p?5] . 

The fundamental character that we can ascribe to two- 
fluids is out of proportion, they are simply dissipative flu- 
ids. Physically: i) The mean free path for both limits are 
very different; ii) Temperature profiles from a transport 
equation, as in the Miiller-Israel-Stewart theory [51]- 
[36] . [37], [35], makes the complete equivalence not vi- 
able at all, enhancing asymmetry. Even so, the master 
fluid idea is attractive, isotropic with free streaming or 
anisotropic with heat flow, to model a number of situa- 
tions; in between we can use the Eddington factor |39| . 
Thus, equivalence in this context can be used as a first 
approximation to study dissipative fluids. 

Our argument in this paper was developed clearly 
thanks to a known approach 27J, to treat matter from a 
comoving frame with the fluid in the local Minkowskian 
frame, called Bondian [30] (usually referred as global or 
Eulerian). We know that the same conclusion can arise 
from Misner-Sharp's comoving coordinates [ID] (local or 
Lagrangian frame) . In this sense, we get an unified treat- 
ment of nonadiabatic fluids. 
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